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INTRODUCTION  John  Thole 

Computation  of  the  shape  of  a  pressurized  thin  spheroid  shell 
composed  of  a  large  number  of  identical  flat  gores  requires  solution  of 
five  nonlinear  partial  differential  equations  of  equilibrium  and  boundary 
conditions  for  three  displacements  of  a  gore.  The  position  of  each  point 
on  the  gore  could  be  determined  as  a  function  of  two  independent  variables, 
namely,  the  distance  of  the  point  from  the  center  line  of  the  gore  and 
the  latitude  of  the  point  with  respect  to  the  equatorial  plane  of  the 
spheroid.  For  an  approximate  solution  to  these  nonlinear  partial  differ¬ 
ential  equations  each  displacement  is  expanded  into  a  power  series  of  the 
distance  of  the  point  from  the  gore  center  line  such  that  the  coefficients 
of  the  power  series  are  at  the  most  functions  of  the  latitude.  By  equating 
to  zero  the  terms  of  the  lowest  power  in  each  equation,  five  new  nonlinear 
differential  equations  for  the  six  coefficients  (two  for  each  displacement) 
are  obtained.  The  approximate  solution  of  these  equations  by  perturbation 
technique  led  to  three  displacement  functions  of  the  gore  with  one  unknown 
constant.  Application  of  variational  method  to  the  total  elastic  strain 
and  potential  energy  of  the  system  led  to  the  evaluation  of  the  above 
mentioned  constant. 

In  this  analysis,  it  is  assumed  that  the  displacement  of  the  spheroid 
due  to  the  gravity  is  negligible  compared  to  the  one  caused  by  the  constant 
internal  pressure.  It  is  also  assumed  that  the  material  of  the  spheroid 
is  homogeneous  and  isotropic  and  it  obeys  Hooke's  law. 
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Finally  this  method  of  analysis  is  applied  to  Echo  II,  the  135 
foot  diameter  passive  communication  satellite.  Since  the  3-ply  aluminum- 
mylar-aluminum  material  of  the  satellite  does  not  satisfy  the  require¬ 
ments  implied  in  the  above  mentioned  assumptions,  therefore  the 
analytical  results  obtained  here  are  expected  to  be  only  an 
approximation  of  the  actual  case. 

THEORY 


In  order  to  determine  the  shape  of  a  thin  spheroid  shell  composed 
of  identical  flat  gores  it  is  sufficient  to  compute  the  true  shape  of 
a  gore  in  the  pressurized  state.  To  do  this,  it  is  assumed  that  the 
material  of  the  spheroid  is  homogeneous  isotropic,  and  linear  elastic. 

It  is  also  assumed  that  the  effect  of  gravity  upon  the  shape  of  the 
spheroid  is  negligible  in  comparison  with  that  of  the  constant  internal 
pressure. 

Let  Cartesian  coordinate  axes  be  drawn  in  such  a  way  that  the 
origin  0  be  located  at  the  center  of  the  balloon,  y  axis  be  along  the 
polar  axis,  and  Z  axis  be  extended  through  the  center  of  a  gore,  as 
shown  in  Figure  1.  Draw  X'Y'Z  axes  from  a  point  P  of  the  gore  in  such 
a  way  that  X'  axis  be  parallel  to  X  axis,  Z'  axis  be  in  the  direction 
OP,  and  Y1  (or  6  )  axis  be  tangent  to  the  meridian  at  point  P  (Figure  1). 

In  the  pressurized  state  of  equilibrium,  the  resultant  of  all  the 
forces  acting  on  an  element  of  the  shell  in  the  direction  X',  0  and  Z' 


Figure  L  -  Coordinate  Axes  and  Displacement  Components 
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vanishes.  Let  N^,  Ng,  and  Nx@  be  nromal  and  shear  tractions  per  unit 
length  of  membrane,  9  be  the  angle  between  the  position  vector  and  the 
equatorial  plane,  and  p  be  the  internal  pressure  of  the  spheroid,  as 


shown  in  Figure  2.  Hence, 


r  10*  +  .  o 

T7 

^  My  6  ^  ^  Nfl  o 

Let  u,  v,  and  w  be  displacement  components  at  a  point  P  in  the  directions 


X',  Y',  and  2'  respectively  (Figure  2).  The  elastic  strains  are  given  by: 
£  ^  -t  L  /  \z 

V  -  -LiV  ^  vvi 


z  -  L  (  .  \  ^  w  \ 

t  ^  ba  a*  r  ; 


Also, 

r N y  =  k  (  +  ft  £&) 

)  si  &  —  k  (  £  q  ^  £x  ) 

(  Ny0~  •<  (  i'"  £*6  (3) 

where  K  =  Et  is  membrane  stiffness,  and  E,  t,  and  u  are  Young's  modulus, 

1-p*  L 

thickness  of  the  shell  and  Poisson's  ratio,  respectively.  Substitution 


of  Eqs .  (2)  into  Eqs.  (3)  leads  to 
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(4) 


Substitution  of  Eqs.  (A)  into  Eqs.  (1)  yields 
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The  boundary  conditions  are: 


'  j\ly  e  -  O  on  y  =  r  ^  Co5  S 

j  ,{J  =*  o/  \A/  CtfS  0  -  o/  V  Six»  £  cw  y  js  r  £t>L s  (6) 

along  the  seam,  where  «*  ^  .  The  first  expression  indicates  that  due 
to  symmetry  shear  stress  vanishes  along  the  seam  where  any  two  adjacent 
gores  meet,  and  the  Second  expression  implies  that  the  points  along  each 
seam  move  in  the  meridional  plane.  Since  the  coefficients  of  only  the 
first  two  terms  of  the  power  series  expansion  of  displacements  will  be 
determined  by  this  method,  continuity  of  slope  in  the  plane  normal  to 
the  seam  cannot  be  used  as  a  third  boundary  condition.  The  use  of  such 


a  boundary  condition  would  lead  to  an  erroneous  solution,  namely,  that 


the  cross  section  of  the  spheroid  with  any  plane  normal  to  its  polar  axis 
is  a  true  circle,  no  matter  what  the  magnitude  of  the  internal  pressure. 
Instead,  a  virtual  displacement  is  given  to  the  pressurized  spheroid  in 
the  state  of  equilibrium  and  change  of  the  total  elastic  strain  and  poten¬ 
tial  energy  of  the  spheroid  and  the  load  is  equated  to  zero. 

Let  r  be  the  radius  of  the  centerline  of  a  gore  before  any  load  is 
applied,  and  a,  b,  c,  d,  f  and  g  be  functions  of  variable  angle  0.  Displace¬ 
ments  W,  U,  and  V  may  be  expanded  into  power  series  as  following: 


\A /  =  <x  r 
U  =  C  y 


(7) 
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Substitution  of  Eqs.  (7)  and  the  last  one  of  Eqs.  (4)  into  Eqs.  (5)  and 


(6)  leads  to 


&A  +  f <)  +  O+/-0  j|  ■+  >z£  ('t'jl>de>+4l:’‘0 
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(8) 

Neglecting  small  terms  of  higher  order  in  Eqs.  (8)  leads  to 
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Since  n,  the  number  of  gores,  is  assumed  to  be  large  therefore  o<  ^  is 
very  small.  Thus  perturbation  technique  can  be  applied  to  the  functions 
a,  b,  c,  d,  f,  and  g  in  the  following  fashion: 

<20  (  6  )  -+■  a,  (0)  -h  '  -  • 

b  (  &  )•=■  b0  (  $  )  of7-  b,  (  & )  ' 

C  (e>)  *  C0  (&)  ■+  ct'C,  (6)  +  - 
dC£>)  =  '<Jo  Ce)  +  ^lcl<  + 
f  Co  )  **  C  & )  -+•  vLf,  (&)  +  — 
i  3-0  J  =  -t-  ^3,  (e)-*-'  ••  (10) 

Substitution  of  Eqs.  (10)  into  Eqs.  (9)  yields  a  set  of  five  nonlinear 
differential  equations  in  terms  of  aQ,  bQ,  c0,  dQ,  fQ,  go,  a^  bx,  c^ 
dl>  etc*  containing  even  powers  of  c<  .  Terms  containing  each  power 
of  o(  in  each  differential  equation  must  independently  vanish.  The 
terms  independent  of  of  must  vanish  in  the  following  way: 

{CAo'+i-pbo'-p  ~  •+■  4-bel-  O 

z1  fe--r^£§  +  ^  =  ° 
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The  solution  of  the  above  equations  could  be  given  as  the  following: 
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2 

The  coefficients  of  oC  in  Eqs.  (9)  must  vanish  in  the  following 
fashion: 

r  ^  1 +  ^ d£-° 

ta,+rc,-^6».,.1BcrYBJi-t 

c !  -  a ,  -+  («+•£)  &]c<ro2e> 

l  a  —  J.  491  —  P-,  dcAi 

ids  &  3a 

The  solution  of  Eqs.  (13)  could  be  written  as: 
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(13) 


(14) 


(15) 


Knowing  the  value  of  Poisson's  ratio,  one  can  easily  evaluate  in  an 
orderly  manner  Y  ,  F-^,  A^,  C-^,  A2,  C2 ,  D^,  and  as  a  function  of 
Constant  B. 


To  evaluate  constant  B,  the  method  of  virtual  displacement  may  be 
utilized.  The  elastic  strain  energy  of  one  quarter  of  a  gore  is  given  by 


T* 


q -  H  Cos© 

v  f 


k  (N*4x  +  Ni  t  -l  rAxAe 


(16) 


Substitution  of  stress  and  strain  values  in  Eq.  (16)  yields 


T-C1j{(^  £  — ^  £  &  -  a-*’ 

7=!  +  £  a|-+ 

±<i§)M  £  ra)+  as- 

■+  h£-  0.  -+  2  \>  —  )Z  i  r  c!  x  d  a 

ro  r  3a  r  -  da'J  <i7) 

In  a  virtual  displacement  corresponding  to  change  of  B  by  j"  B  the  work 
done  by  the  internal  pressure  for  a  quarter  of  a  gore  is  given  by 


^  W 


£e>«  t 
1 


p 

f  xt 

1  J 

J 
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in  which  ft  is  the  total  work  done  by  the  internal  pressure  for  a  quarter 


(18) 


gore,  and  W  is  the  radial  displacement  given  by  Eqs.  (7).  Hence,  the 
Constant  B  could  be  determined  from 

r  £  r  x=  r*  9 


It  _  . 


6=0 


x=  0 


^rAxd  6 


(19) 


Thus  for  a  given  pressure  p,  Eq.  (19)  in  conjunction  with  Eq.  (17),  (7), 


(10),  (12),  (14)  and  (15)  furnishes  B  which  upon  substitution  into  Eq. 
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(15)  gives  coefficients  A^,  A£,  etc.  Substitution  of  these  coefficients 
into  Eq.  (14)  in  conjunction  with  Eqs.  (7),  (10)  and  (12)  gives  the 
three  displacement  functions  u,  v,  and  w. 

APPLICATION  OF  THEORY  TO  ECHO  II 

The  135  fool:  diameter  Echo  II  is  composed  of  106  identical  gores 
with  a  maximum  width  of  4  feet  which  occurs  in  the  equatorial  plane  of 
the  balloon.  The  3-ply  material  of  the  balloon  is  composed  of  0.00035 
in.  thick  mylar  sandwiched  between  two  layers  of  aluminum  of  0.00018  in. 
thickness  each.  The  adjoining  gores  are  located  edge  to  edge  with  each 
other  with  one  inch  wide  tape  of  the  same  material  sealing  the  seams 
lengthwise.  When  separation  of  the  two  half -canisters  takes  place  in 
the  initial  phase  cf  orbiting,  the  residual  air  and  water  vapor  inside  the  folded 
balloon  which  is  of  the  order  of  magnitude  of  one  millimeter  inflates  the 
satellite.  While  initially  flat,  each  gore  forms  a  narrow  transverse 
portion  of  a  half  circular  cylindrical  surface  of  135  feet  diameter. 

Thus,  the  initial  shape  of  the  balloon  fc  a  spheroid  with  0.356  in.  maximum 
deviation  from  a  sphere  which  occurs  in  the  equatorial  plane. 

Due  to  the  solar  radiation  skin  temperature  increases  and  as  a 
result  inflation  material  sublimes  and  builds  up  pressure  until  a 
maximum  internal  pressure  of  225  microns  is  reached.  This  maximum 
pressure  corresponds  to  a  membrane  stress  of  1.7  lb/in. 

Although  the  material  of  the  balloon  is  rather  nonhomogeneous  and 
anisotropic  and  it  does  not  quite  obey  Hooke's  law,  nevertheless  for  the 
purpose  of  approximation,  it  is  assumed  that  the  requirements  of  homogeneity, 

isotropy  and  linear  elasticity  necessitated  by  this  analysis  are  met. 
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The  average  value  of  K  corresponding  to  1.70  lb/in.  skin  traction  (given 
by  Tests  No.  1,  3,  4,  and  6  of  Reference  3)  is  1050  lb/in. 

Since  nominal  Young's  modulus  for  aluminum  is  10.5  x  10^  psi  and 
for  mylar  is  .55  x  10^  psi,  therefore  it  is  expected  that  aluminum 
should  carry  most  of  the  load.  The  nominal  value  of  K  for  aluminum 
layers  is  given  by: 


Kn 


Et 


r 


4240  lb/in. 


in  which  ^  =  .33  has  been  inserted.  This  value  of  Kn  is  much  larger 
than  the  experimental  value  K  *  1050  lb/in.  Therefore,  the  actual  Young's 
modulus  and  perhaps  Poisson's  ratio  for  thin  layers  of  aluminum  are 
smaller  than  the  nominal  ones. 

Following  the  steps  suggested  in  the  last  paragraph  of  the  previous 

section  one  could  get  an  Eigen  value  equation  for  B  whose  coefficients 

are  functions  of  parameter  m  -  Prg  •  For  different  values  of  maximum 

K  o( 

pressure  the  fundamental  Eigen  value  for  B  is  computed  and  listed  in 

Table  1.  Hence  for  p  *  225  micron  (or  .000435  lb/in^). 

.000435  x  810 


m 


1050  x 


'106  1 


3.81 


which  corresponds  to  B  =  -.319  in  Table  1.  Thus,  the  maximum  differential 
radial  displacement  at  the  midgore  and  at  the  seam  which  occurs  in  the 
equatorial  plane  is  given  by 


\A/ 


-  w 

1  -  13 

x2 

b~o 

O 

» 

r 

>=  0 

X.2V 

1 

=  O.ZZ7 


This  corresponds  to  a  maximum  radial  difference 

A  =  0.356  -  0.227  =  0.129  in  which  0.356  is  in  its  initial 


unpressurized  value. 
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TABLE  I 


Values  of  B  versus  m 


m 

B 

m 

B 

m 

B 

m 

B 

m 

B 

0 

0 

2.0 

-.294 

4.0 

-.320 

6.0 

-.320 

8.0 

-.317 

.1 

-.043 

2.1 

-.298 

4.1 

-.320 

6.1 

-.320 

8.1 

-.317 

.2 

-.079 

2.2 

-.300 

4.2 

-.320 

6.2 

-.320 

8.2 

-.317 

.3 

-.109 

2.3 

-.303 

4.3 

-.320 

6.3 

-.319 

8.3 

-.317 

.4 

-.135 

2.4 

-.305 

4.4 

-.320 

6.4 

-.319 

8.4 

-.317 

.5 

-.157 

2.5 

-.307 

4.5 

-.320 

6.5 

-.319 

8.5 

-.317 

.6 

-.176 

2.6 

-.309 

4.6 

-.320 

6.6 

-.319 

8.6 

-.316 

.7 

-.193 

2.7 

-.311 

4.7 

-.320 

6.7 

-.319 

8.7 

-.316 

.8 

-.207 

2.8 

-.312 

4.8 

-.321 

6.8 

-.319 

8.8 

-.316 

.9 

-.220 

2.9 

-.313 

4.9 

-.321 

6.9 

-.319 

8.9 

-.316 

1.0 

-.213 

3.0 

-.314 

5.0 

-.321 

7.0 

-.319 

9.0 

-.316 

1.1 

-.241 

3.1 

-.315 

5.1 

-.321 

7.1 

-.318 

9.1 

-.316 

1.2 

-.250 

3.2 

-.316 

5.2 

-.321 

7.2 

-.318 

9.2 

-.316 

1.3 

-.258 

3.3 

-.317 

5.3 

-.321 

7.3 

-.318 

9.3 

-.316 

1.4 

-.265 

3.4 

-.318 

5.4 

-.320 

7.4 

-.318 

9.4 

-.316 

1.5 

-.272 

3.5 

-.318 

5.5 

-.320 

7.5 

-.318 

9.5 

-.315 

1.6 

-.277 

3.6 

-.319 

5.6 

-.320 

7.6 

-.318 

9.6 

-.315 

1.7 

-.282 

3.7 

-.319 

5.7 

-.320 

7.7 

-.318 

9.7 

-.315 

1.8 

-.287 

3.8 

-.319 

5.8 

-.320 

7.8 

-.317 

9.8 

-.315 

1.9 

-.291 

3.9 

-.320 

5.9 

-.320 

7.9 

-.317 

9.9 

-.315 
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CONCLUSIONS 


The  maximum  differential  radial  displacement  of  a  pressurized 
spheroid,  which  is  composed  of  identical  flat  gores,  in  a  zero  gravity 
field,  computed  by  this  technique  is  believed  to  be  a  very  good  approxi 
mation  to  the  true  one,  provided  that  the  material  is  homogeneous,  iso¬ 
tropic,  and  linearly  elastic. 
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